Combining the modified discrete element method with the virtual element method for fracturing of porous media     

Halvor Møll Nilsen  Idar Larsen  Xavier Raynaud 《Computational Geosciences》2017,21(5-6):1059-1073

Simulation of fracturing processes in porous rocks can be divided into two main branches: (i) modeling the rock as a continuum enhanced with special features to account for fractures or (ii) modeling the rock by a discrete (or discontinuous) approach that describes the material directly as a collection of separate blocks or particles, e.g., as in the discrete element method (DEM). In the modified discrete element (MDEM) method, the effective forces between virtual particles are modified so that they reproduce the discretization of a first-order finite element method (FEM) for linear elasticity. This provides an expression of the virtual forces in terms of general Hook’s macro-parameters. Previously, MDEM has been formulated through an analogy with linear elements for FEM. We show the connection between MDEM and the virtual element method (VEM), which is a generalization of FEM to polyhedral grids. Unlike standard FEM, which computes strain-states in a reference space, MDEM and VEM compute stress-states directly in real space. This connection leads us to a new derivation of the MDEM method. Moreover, it enables a direct coupling between (M)DEM and domains modeled by a grid made of polyhedral cells. Thus, this approach makes it possible to combine fine-scale (M)DEM behavior near the fracturing region with linear elasticity on complex reservoir grids in the far-field region without regridding. To demonstrate the simulation of hydraulic fracturing, the coupled (M)DEM-VEM method is implemented using the Matlab Reservoir Simulation Toolbox (MRST) and linked to an industry-standard reservoir simulator. Similar approaches have been presented previously using standard FEM, but due to the similarities in the approaches of VEM and MDEM, our work provides a more uniform approach and extends these previous works to general polyhedral grids for the non-fracturing domain.  相似文献   

Architecture and early evolution of the Oslo Rift     

B. Sundvoll  B.T. Larsen   《Tectonophysics》1994,240(1-4):173-189

A revised assessment of architecture and pre-rift fabric connections of the Oslo Rift has been undertaken and linked to a new appraisal of observations and data related to the initial phase of the rift evolution. In addition to half-graben segmentation, accommodation zones and transfer faults are readily identified in the linking sectors between the two main grabens and between graben segments. Axial flexures are proposed between facing half-grabens. The accommodation zones were generally sites of volcanism during rifting. Pre-rift tectonic structures played an influential role in the rift location and development. The deviant N-S axis of the Vestfold graben segment is viewed as related to pre-rift structural control through faults and shear zones. This area was probably a site of Proterozoic/Palaeozoic crustal and lithospheric attenuation.

Field evidence suggests that the rift started as a crustal sag with no apparent surface faulting in a flat and low-lying land at a time about 305–310 Ma. Volcanism, sub-surface sill intrusion and faulting started about simultaneously some time after the initial sag (300–305 Ma). Faulting and basaltic volcanism were initially localized to transfer faults along accommodation zones and a NNW-SSE transtensional zone along the eastern margin of the incipient Vestfold graben segment. This transtensional zone was probably created by right-lateral simple shear tracing pre-rift structures in response to a regional stress field with the tensional axis normal and the maximum compressional axis parallel to the NNE-SSW-trending rift axis.  相似文献   

The role of soil processes in determining mechanisms of slope failure and hillslope development in a humid-tropical forest eastern Puerto Rico   总被引：1，自引：0，他引：1  

Andrew Simon  Matthew C. Larsen  Cliff R. Hupp 《Geomorphology》1990,3(3-4)

Translational failures, with associated downslope earthflow components and shallow slides, appear to be the primary mechanism of hillslope denudation in the humid tropical forests of the mountains of eastern Puerto Rico. In-situ weathering of quartz diorite and marine-deposited volcaniclastics produces residual soil (saprolite; up to 21 m deep) / weathered rock profiles. Discontinuous zones of contrasting density and permeability particularly in quartz-diorite slopes at 0.5 m, and between 3 and 7 m, create both pathways and impedances for water that can result in excess pore pressures and, ultimately, aid in determining the location of failure planes and magnitudes of slope failures. In combination with relict fractures which create planes of weakness within the saprolite, and the potential significance of tensile stresses in the upper zone of saprolite (hypothesized to be caused by subsurface soil creep), shear failure can then occur during or after periods of heavy rainfall.Results of in-situ shear-strength testing show negative y-intercepts on the derived Mohr-Coulomb failure envelopes (approximately 50% of all tests) that are interpreted as apparent tensile stresses. Observation of tension cracks 1–2 m deep support the test data. Subsurface soil creep can cause extension of the soil and the development of tensile stresses along upper-slope segments. Shear-strength data support this hypothesis for both geologic types. Apparent values of maximum and mean tensile stress are greatest along upper slopes (16.5 and 6.29 kPa). Previously documented maximum rates of downslope movement coincided with local minima of shear strength, and the shear-strength minimum for all tests was located near 0.5 m below land surface, the shallow zone of contrasting permeabilities. These results indicate that subsurface soil creep, a slow semi-continuous process, may exert a profound influence on rapid, shallow slope failures in saprolitic soils.Data indicate that cove slopes in quartz diorite tend to be the most unstable when saturation levels reach 75%. Deep failures (7 m deep) appear the most critical but not the most frequent because pore pressure build-up will occur more rapidly in the upper perched zone of translocated clays before reaching the lower zone between 3 and 7 m. Frequent shallow failures could reduce the probability of deeper failures by removing overburden and reducing shear stress at depth. Deep failures are more likely to result from storm events of great duration and intensity.Sixty-six ‘naturally occurring’ and more than 100 ‘road-related’ landslides were mapped. Forest elevations exceed 1000 m, but the majority of these failures were found between 600 and 800 m in elevation. This appears to be the area where there is sufficient concentration of subsurface water to result in excess pore pressures. The high percentage of slope failures in the 600–800-m range, relative to the percentage at higher elevations, suggests that differences in soil-water processes are responsible for the form of these mountain slopes. Steep linear segments are maintained at higher elevations. Slope angles are reduced in the 600–800-m range by frequent shallow slides, creating a largely concave surface. In combination, slope segments above 800 m, and those between 600 and 800 m, produce the characteristic form of the mountains of eastern Puerto Rico.  相似文献   

Origin of high manganese concentrations in coal mine drainage, eastern Tennessee     

Daniel Larsen  Richard Mann 《Journal of Geochemical Exploration》2005,86(3):143-163

The origin of high dissolved manganese concentrations in slightly acidic mine runoff from a surface mine operated by the Cumberland Coal Company in eastern Tennessee was investigated. Mineralogical and chemical analyses were performed on 31 samples of sandstone, shale, coal, and mudstone from the mine to identify the sources and stratigraphic distribution of high extractable manganese contents in the spoil materials. The samples were analyzed for their bulk mineral content by X-ray diffraction, net acid-base potential, and reaction to 2 or 4 chemical extraction procedures. A limited number of samples were analyzed for petrographic characteristics, clay mineral composition by X-ray diffraction, and mineral compositions by electron microprobe. Analysis of the data and consideration of the geochemical conditions at the mine were used to identify probable sources for the high extractable manganese contents.The results indicate 2 prominent, independent sources of extractable manganese. The first source is exchangeable manganese on clay minerals (mainly illite + muscovite and chlorite) and is concentrated in shale and mudstone rock types. The second and more significant source is manganese in siderite concretions and cement, mainly in shale and mudstone. Comparison to other coal-bearing strata indicates that manganese-rich siderite is common in fresh- to brackish-water subaqueous sediments that overlie coal. This is especially the case for coals formed in wet, tropical environments.Ratios of manganese to calcium and magnesium in mine runoff suggest that manganese from siderite is the major cause of the high dissolved manganese contents. A conceptual model is developed to explain the high manganese contents of the mine runoff. Oxidation of pyrite creates mildly acidic waters that are subsequently partially neutralized by reaction with impure siderite. Solubilized manganese remains dissolved in the slightly acidic runoff water, whereas dissolved iron precipitates as ferric hydroxide or goethite. Consideration of data from other coal mining regions suggests that similar reactions involving impure siderite may be responsible for high manganese concentrations in acidic to slightly acidic mine runoff. Geochemical reaction path modeling of pyrite and impure siderite with rainwater illustrate how resulting water compositions may vary depending on pyrite to siderite ratios in spoil materials. Spoil water compositions from the Cumberland mine are largely consistent with reaction of pyrite and impure siderite in proportions observed in the sediments; however, deviations may be explained by minor mixing with waters that reacted only with impure siderite or clay mineral exchange reactions.  相似文献   

Gas-Phase OH Oxidation of Monoterpenes: Gaseous and Particulate Products     

Bo. R. Larsen  Dario Di Bella  Marianne Glasius  Richard Winterhalter  Niels R. Jensen  Jens Hjorth 《Journal of Atmospheric Chemistry》2001,38(3):231-276

Smog chamber experiments have beenconducted in which cyclic monoterpenes were oxidisedin the gas phase by OH. The evolved secondary organicaerosol (SOA) was analysed by LC-MSⁿ and thegas-phase products were analysed by FT-IR. Theconcentrations of the identified compoundscorresponded to carbon mass balances in the range of40%–90%. The identified compounds in the particularphase corresponded to 0.5%–4.2% of the reactedcarbon. The most abundant compounds in SOA fromterpenes with an endocyclic C=C double bond wereC₁₀-keto-aldehydes, C₁₀-keto-carboxylicacids, C₁₀-hydroxy-keto-carboxylic acids, andC₁₀-hydroxy-keto-aldehydes (pinonaldehyde,pinonic acid, hydroxy-pinonic acid isomers, andhydroxy-pinonaldehyde isomers from -pinene;3-caronaldehyde, 3-caronic acid, hydroxy-3-caronicacid isomers, and hydroxy-3-caronaldehyde isomers from3-carene). The most abundant compounds in SOA fromterpenes with an exocyclic C=C double bond wereC₉-ketones, C₉-dicarboxylic acids, andC₁₀-hydroxy-keto-carboxylic acids (nopinone,pinic acid, and hydroxy-pinonic acid isomers from-pinene; sabinaketone, sabinic acid andhydroxy-sabinonic acid isomers from sabinene).Decarboxylated analogues of most of the compounds werepresent in SOA in minor concentrations, such asC₉-keto-carboxylic acids (norpinonic acid,nor-3-caronic acid) and C₈-dicarboxylic acids(norpinic acid, nor-3-caric acid, norsabinic acid). InSOA from limonene, which contains an endocyclic aswell as an exocyclic C=C double bond, the mostabundant compounds were a C₁₀-keto-aldehyde andits oxo-derivative (limononaldehyde and keto-limononaldehyde) together with hydroxy-derivatives of aC₁₀-keto-carboxylic acid (isomers ofhydroxy-limononic acid). Also aC₁₀-keto-carboxylic acid (limononic acid) waspresent together with minor concentrations of aC₉-dicarboxylic acids (limonic acid), itsoxo-derivative (keto-limonic acid), and itsdecarboxylated analogue (norlimonic acid). Mechanisticpathways for the formation of these products, some ofwhich are identified here for the first time, areproposed.  相似文献   

The last glaciation of the Arctic volcanic island Jan Mayen     

Astrid Lyså  Eiliv A. Larsen  Johanna Anjar  Naki Akçar  Morgan Ganerød  Asbjørn Hiksdal  Roelant Van Der Lelij  Christof Vockenhuber 《Boreas: An International Journal of Quaternary Research》2021,50(1):6-28

The volcanic island of Jan Mayen, remotely located in the Norwegian-Greenland Sea, was covered by a contiguous ice cap during the Late Weichselian. Until now, it has been disputed whether parts of the island south of the presently glaciated Mount Beerenberg area were ever glaciated. Based on extensive field mapping we demonstrate that an ice cap covered all land areas and likely also extended onto the shallow shelf areas southeast and east of the island. Chronological interpretations are based on K-Ar and ⁴⁰Ar/³⁹Ar dating of volcanic rocks, cosmogenic nuclide (³⁶Cl) surface exposure dating of bedrock and glacial erratics, and radiocarbon dating. We argue that ice growth started after 34 ka and that an initial deglaciation started some 21.5–19.5 ka in the southern and middle parts of the island. In the northern parts, closer to the present glaciers, the deglaciation might have started later, as evidenced by the establishment of vegetation 17–16 cal. ka BP. During full glaciation, the ice cap was likely thickest over the southern part of the island. This may explain a seemingly delayed deglaciation compared with the northern parts despite earlier initial deglaciation. In a broader context, the new knowledge of the Late Weichselian of the island contributes to the understanding of glaciations surrounding the North Atlantic and its climate history.  相似文献   

Geochemical and isotopic evidence for upward flow of saline fluid to the Mississippi River Valley alluvial aquifer,southeastern Arkansas,USA     

Larsen  Daniel  Paul  Justin  Cox  Randy 《Hydrogeology Journal》2021,29(4):1421-1444

Hydrogeology Journal - Groundwater from the Quaternary Mississippi River Valley Alluvial (MRVA) aquifer in southeastern Arkansas (SE AR), USA, has higher salinity compared to other MRVA...  相似文献   

Origin and extent of fresh groundwater, salty paleowaters and recent saltwater intrusions in Red River flood plain aquifers, Vietnam   总被引：1，自引：0，他引：1  

Luu T. Tran  Flemming Larsen  Nhan Q. Pham  Anders V. Christiansen  Nghi Tran  Hung V. Vu  Long V. Tran  Hoan V. Hoang  Klaus Hinsby 《Hydrogeology Journal》2012,20(7):1295-1313

A model has been established on the origin and extent of fresh groundwater, salty paleowaters and saltwater from recent seawater intrusions in the Red River flood plain in Vietnam. This was done with geological observations, geophysical borehole logging and transient electromagnetic methods. Salt paleowater is present up to 50–75?km from the coastline, with occurrence controlled by the Holocene transgression. A density-driven leaching of salty porewater has occurred from high-permeability Holocene sediments into underlying Pleistocene deposits, whereas diffusion has dominated in low-permeability layers. In the Pleistocene aquifer, the highest content of dissolved solids is found below two intrinsic valleys with Holocene marine sediments and along the coastline. Recent intrusion of saltwater from the South China Sea is observed in shallow groundwater 35?km inland, probably a result of transport of salty water inland in rivers or leaching of paleowaters from very young near-coast marine sediments. The observed inverted salinity profile, with high saline water overlying fresher groundwater, has been formed due to the global eustatic sea-level changes during the last 8,000–9,000?years. The proposed model may therefore be applicable to other coastal aquifers, with a proper incorporation of the local geological environments.  相似文献   

A 10Be‐based reconstruction of the last deglaciation in southern Sweden     

Johanna Anjar  Nicolaj K. Larsen  Lena Håkansson  Per Möller  Henriette Linge  Derek Fabel  Sheng Xu 《Boreas: An International Journal of Quaternary Research》2014,43(1):132-148

We present 23 cosmogenic surface exposure ages from 10 localities in southern Sweden. The new ¹⁰Be ages allow a direct correlation between the east and west coasts of southern Sweden, based on the same dating technique, and provide new information about the deglaciation of the Fennoscandian Ice Sheet in the circum‐Baltic area. In western Skåne, southernmost Sweden, a single cosmogenic surface exposure sample gave an age of 16.8±1.0 ka, whereas two samples from the central part of Skåne gave ages of 17.0±0.9 and 14.1±0.8 ka. Further northeast, in southern Småland, two localities gave ages ranging from 15.2±0.8 to 16.9±0.9 ka (n=5) indicating a somewhat earlier deglaciation of the area than has previously been suggested. Our third locality, in S Småland, gave ages ranging from 10.2±0.5 to 18.4±1.6 ka (n=3), which are probably not representative of the timing of deglaciation. In central Småland one locality was dated to 14.5±0.8 ka (n=3), whereas our northernmost locality, situated in northern Småland, was dated to 13.8±0.8 ka (n=3). Samples from the island of Gotland suggest deglaciation before 13 ka ago. We combined the new ¹⁰Be ages with previously published deglaciation ages to constrain the deglaciation chronology of southern Sweden. The combined deglaciation chronology suggests a rather steady deglaciation in southern Sweden starting at c. 17.9 cal. ka BP in NW Skåne and reaching northern Småland, ～200 km further north, c. 13.8 ka ago. Overall the new deglaciation ages agree reasonably well with existing deglaciation chronologies, but suggest a somewhat earlier deglaciation in Småland.  相似文献   

Late Quaternary biotic and abiotic controls on long-term sediment flux in a northern Australian tropical river system     

Annegret Larsen  Jan-Hendrik May  Xavier Carah 《地球表面变化过程与地形》2019,44(12):2494-2509

The modern distribution of monsoonal rainforest in the Australian tropics is patchy and is mainly associated with river corridors and groundwater springs, which indicates a strong dependence on hydrologic and geomorphic conditions. While their present distribution is well known, very little data exists on past spatial and temporal dynamics of these ecosystems, or their medium- to longer-term controls. Factors such as (i) fire frequency and type, and/or (ii) hydroclimatic conditions (e.g. droughts) have been proposed to control riverine corridor rainforest extent. Recent observations, however, also suggest an additional (iii) geomorphic control induced by alluvial knickpoint migration. Sediment sequences provide valuable archives for the reconstruction of longer-term (a) floodplain sedimentary dynamics, (b) local vegetation history, and (c) catchment-wide fire histories. This study investigates such a sediment sequence at Wangi Creek, and shows that a phase of aggradation, lasting ~4000 years, was recently disrupted by channel incision and floodplain erosion. The aggradational phase is characterized by sand deposition with average vertical floodplain accretion rates of 0.8 cm/yr and includes phases of soil development. The recent incisional phase has changed hydro-geomorphic conditions and caused widespread degradation of vegetation, erosion and lowering of the macro-channel surface. While there is no evidence in our data for an erosional event of similar magnitude since the onset of late Holocene floodplain aggradation, Wangi Creek experienced significant erosion and incision immediately before ~4000 years, providing the first evidence for a tropical cut-and-fill river system. We hence argue that phases of aggradation mainly controlled by biotic processes alternate and depend on feedbacks with incision phases controlled mainly by abiotic processes. The results show that eco-hydro-geomorphic feedbacks may play a crucial role in the medium- to longer-term history of tropical fluvial systems and need to be considered when interpreting fluvial archives with regards to climate, fire or human induced change. © 2019 John Wiley & Sons, Ltd.  相似文献